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Active metamaterials: sign of refraction index and gain-assisted dispersion 

management 
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We derive an approach to define the causal direction of the wavevector of modes in optical meta- 
materials, which in turn, determines signs of refractive index and impedance as a function of real 
and imaginary parts of dielectric permittivity and magnetic permeability. We use the developed 
technique to demonstrate that the interplay between resonant response of constituents of metama- 
terials can be used to achieve efficient dispersion management. Finally we demonstrate broadband 
dispersion-less index and impedance matching in active nanowire-based negative index materials. 
Our work opens new practical applications of negative index composites for broadband lensing, 
imaging, and pulse-routing. 



Recent research in the area of negative index mate- 
rials (NIMs)i has resulted in a number of exciting ap- 
plications, including super-imaging^ and subwavelength 
light compression 4 . However, the majority of these appli- 
cations suffer from substantial material absorption 5 - and 
frequency dispersion of NIM composites. Optical gain 
has been suggested to minimize, and potentially elimi- 
nate absorption losses 6 -^. While effect of gain on prop- 
agation length of optical signals is straightforward, its 
effect on material dispersion has not been completely 
understood. Furthermore, the very question of sign 
of refractive index in active metamaterials is somewhat 
controversia l 10 ' 15 . In this Letter we present a univer- 
sal approach for imposing causality in active and passive 
materials, and use this technique to analyze the perspec- 
tives of gain-assisted dispersion management beyond loss 
compensation. 

Our results can be used to determine the sign of re- 
fractive index of active or passive media, as well as 
in a number of analytical and/or numerical solutions 
of Maxwell equations relying on plane-wave representa- 
tions. Our technique is illustrated using an example of 
nanowire-based NIM structure, originally proposed inii 
and experimentally realized in 12 . It is shown that fre- 
quency independent (negative) index and impedance can 
be achieved in the same metamaterial with position- 
dependent gain in weak gain regime. A combination 
of broadband impedance and refractive index has a po- 
tential to open new exciting applications of dispersion- 
managed NIMs, in broadband optical processing, packet 
routing, and non-reflective lensing. 

Index of refraction n p h is one of the most fundamen- 
tal optical properties of the media. Its magnitude relates 
the magnitude of wavevector k of plane electromagnetic 
wave to frequency u>: |n p h| = \k\c/u), and thus describes 
the phase velocities of waves in the material—. In par- 
ticular, riph enters the equations for reflectivity, Dopplcr 
effect, and some nonlinear phenomena. Apart from n p h, 
reflectivity of a material also depends on its impedance 
Z . For an isotropic material with (complex) dielectric 



permittivity e and magnetic permeability fx, n p h and Z 
are calculated via 1 ^: 



n ph = ±^/eJl 
Z = n/riph. = ±\/Ji/e, 



(1) 
(2) 



Note that while the magnitude of n p h and Z are 
completely determined by material paramaters (e 
and /tt)i£, their signs have recently instigated some 
controvers y 9 ' 10 ' 14 ' 15 , which can be traced to different 
treatments of causality principle. Moreover, the authors 
of£ suggest that Maxwell equations can be solved cor- 
rectly regardless the selection of sign of refractive index. 
Such a freedom of choice, however is accompanied by 
a requirement to adjust the signs in equations describ- 
ing phase velocity- related phenomena, e.g. Snclls law 9 , 
and still require imposing causality (identical to^) when 
solving Maxwell equations. 

From a mathematical standpoint, imposing causal- 
ity principle is equivalent to selecting the sign of the 
wavevector of a plane wave propagating away from the 
source. Here we assume that such a propagation takes 
place along the positive z direction and therefore focus 
on the k z component of the wavevector. The authors 
O g ,io 1 i4 p r0 p 0se t select the sign of k z , enforcing posi- 
tive direction of Poynting vector (associated with energy 
flux). The authors o£i£ suggest that causality requires 
exponential decay (k" > Qi£) of waves propagating in- 
side passive materials and exponential growth [k" < 0) 
of waves propagating inside active media. 

While all causality requirements discussed 
abov o 9 ' 10 ' 14 ' 15 coincide for the case of passive ma- 
terials, they are not directly applicable for active media 
and are therefore not universal. Indeed, enforcing the 
sign of energy flux is physical only in transparent mate- 
rials. Materials with opposite signs of e and fi (known 
as single-negative materials^ 3 .) reflect the majority 
of incident radiation. Enforcing decay/growth of field 
based solely on passive/active state of material yields 
nonphysical results (such as abrupt disappearance of 
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surface plasmon polariton waves) when media undergoes 
smooth transition from low-loss to low-gain stated. 

To resolve the above controvers y 10 ^ 5 , we propose to 
simultaneously consider "eigen transparency" of the ma- 
terial (its transparency in the absence of losses and gain, 
e" = fi" = 0) along with absorption (or gain) state of 
the material. Clearly, electromagnetic radiation should 
decay inside all passive medial. It should also grow in- 
side transparent (double-negative or double-positive) ac- 
tive materials. Non-transparent (single-negative) mate- 
rials do not support propagating modes, and thus should 
reflect all incident light. Energy can penetrate these 
structures only in the form of exponentially decaying 
wavesii. Since decay/growth of EM waves can be re- 
lated to the sign of imaginary part of refraction index, 
our arguments, summarized in Table I, provide complete 
solution to the problem of selection of direction of the 
wavevector of plane waves. For isotropic media, the de- 
veloped technique also provides a solution to selection of 
the sign of n' ph , which should be identical to that of k' z , 
yielding "conventional" Snell's lawi 

For passive media, our results agree with those of^^, 
and withi 5 -, relying on the pre-selected branch cut in the 
complex plane when calculating the square root in Eq.([T]). 
We note however, that Table I cannot be reduced to such 
a cut. Indeed, an optical material can fall into one of the 
four cases: it has either negative (n' ph < 0)or positive 
(n' ph > 0) refractive index, and it attenuates (n ph > 0) 
or amplifies {n ph < 0) incoming radiation^. Selection of 
any single complex plane cut in Eq. ([1]) would immedi- 
ately limit the number of possible {^phJ n phi combina- 
tions to two, and therefore in general is not correct. 

While it is impossible to determine a universal single 
cut in the complex plane for n p h, requirements of Table 
I can be formally satisfied by the following procedure: 
starting from material parameters e and /i, one first cal- 
culates y'e and *//!, cutting the complex plane along neg- 
ative imaginary axis (see^ and discussion in Ref.— on 
implications of different selection of signs of k" z in non- 
magnetic media). Refraction index and impedance are 
then calculated as n p h = \fe ■ y/ji) Z = ^/Ji/\/e 

We note that the above procedure can be general- 
ized to other classes of materials and excitation waves. 
For example, similar transparency/active state argu- 
ments can be employed to find the "handedness" (re- 
lationship between directions of E, H, and k) of modes 
in active anisotropy-based waveguides^, as well as the 
handedness of waves in active media excited by evanes- 
cent radiatio n 16 ! 20 (the latter case is typically realized in 
gain-assisted compensation of losses of surface plasmon 
polaritonsi). 

We now employ the developed technique to analyze the 
gain-assisted dispersion management in active negative- 
index metamaterials. To illustrate our approach, we 
select nanowire- based optical NIM syste m 11 ! 12 . The 
meta-atom of this composite comprises two plasmonic 
nanowires, schematically shown in inset of FigQJ). As 
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TABLE I: The direction of wavevector (and thus the sign of 
refractive index) in optical material is related to the interplay 
between transparency and gain/loss state of the media. The 
table summarizes this dependence. First column represents 
the transparency state of the material, determined by the 
sign of product e'fj/; The sign of \e\e" + \fi\/J." (second column) 
determines whether material is passive ( |e|e" + > 0) or 

active ( |e|e" + Mm" < 0)- The sign of the refractive index 
is selected to satisfy the requirement for wave attenuation 
(fc" > 0) or growth (fc" < 0) (third column) 



described in detail inii, dielectric and magnetic response 
of such a metamaterial can be related to symmetric and 
anti-symmetric combinations of surface plasmon polari- 
tons excited in nanowire pairs. In the limit of small con- 
centration of nanowires p, effective dielectric permittivity 
and magnetic permeability of such a mix can be qualita- 
tively described byii: 



e c ff = £h 
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(3) 



(4) 



where, r, I, and d correspond to nanowire ra- 
dius, length, and separation between two wires, 
and remaining parameters are given by: fl 2 — 
k 2 l 2 (\n[l/2r]+i^kl/2)/(A In [1 + e h l/2r]), C 2 = 
e h /(41n[d/r]), g 2 - k 2 e h [l + */(2A 2 f[A] ln[d/r])] , A = 
W 3 ^, f(A) = (l-i)Ji[(l + i)A]/(AJ [(l + i)A]), 




A, urn 



FIG. 1: Effective permittivity (a), permeability (b), refractive 
index (c) and impedance (d) of the passive nanowire NIMs 
described in text and schematically shown in inset in (b); 
solid lines corresponds to A = 0.04, dashed to A = 0, dash- 
dotted to A — 0.0175. Inset in (a) shows real and imaginary 
parts of permittivity of host material 
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with k = 2n/X = lu/c, A is being wavelength in the vac- 
uum, and e m and e/, being permittivities of nanowires 
and host materials. Here we assume that e m of silver 
wires is described by the Drude model3i, and further as- 
sume that host dielectric consists of a polymer (eo — 1.5) 
doped with quantum dots^ 2 -, qualitatively described by 
Lorentz model: 



where loq is the resonant frequency 7 is the damping con- 
stant, and A is the macroscopic analog of Lorentz oscilla- 
tor strength, which formally describes gain in the system 
and can be related to the concentration of quantum dots 
and the fraction of quantum dots in excited state. A > 
corresponds to lossy materials; A = represents the case 
when the number of excited quantum dots is equal to the 
number of dots in the ground state; A < corresponds 
to the inverted (gain) regime™. The permittivity of the 
host medium and corresponding permittivity of the NIM 
system for different pump rates are shown in inset of 
Figffli. 

Figdk, illustrates perspectives of dispersion manage- 
ment in lossy (A > 0) nanowire composites with p = 0.1, 
r = 25nm, I = 700nm, d — 120nm, and Ao = 2ttc/ujq = 
1.5/im 7 = 0.628/im -1 . It is clearly seen that disper- 
sion of host media can completely compensate the dis- 
persion of refractive index and impedance of the NIM 
system. Note, however, that due to their different depen- 
dence on e and /1, broadband refractive index and broad- 
band impedance are realized at different values of oscil- 
lator strength A in the single-oscillator model assumed 
here. We suggest that benefits of impedance-matching 
can be combined with benefits of index-matching in the 
same system where A is (adiabatically) changed from 
A ~ 0.0175 corresponding to dZ/duu — at the inter- 
face to A ~ 0.04 corresponding to dn-ph/du) = in the 
core of the system. In quantum dot materials, spatial 
change of A can be easily achieved by changing quantum 
dot doping or external pumping rate. 

Although passive host (A > 0) does yield broadband 
frequency-independent n p h and Z, it also increases total 
absorption in NIM structure, further limiting its practical 
size to < 1 • • -10/im 5 . Active quantum dots, on the other 
hand can simultaneously reduce absorption in the sys- 
tem and provide versatile gain-assisted dispersion man- 
agement. Note, that such a modulation of n p h or Z does 
not require full compensation of propagation losses. 

Gain-assisted dispersion management in active 
nanowire composites with A < 0, Ao = lA/im, 
7 = 0.129//m _1 , p = 0.09, r = 25nm, I = 720nm, and 
d = 120nm are shown in FigfJ] Note that refraction 
index, impedance, as well as reflectivity between vac- 
uum and nanowire metamaterial are continuous when 
material switches between active and passive states. In 
contrast to this behavior, transition between transparent 
and "metallic" regimes yields a discontinuity in reflectiv- 
ity (the discontinuity disappears when thickness of gain 




FIG. 2: Real and imaginary parts of the effective permittivity 
(a,b), permeability (e,f), refractive index (c,d), real part of 
impedance (g) of of active nanowire NIMs, and reflection from 
the semi-infinite slab of this system (i). Solid, dashed, and 
dash-dotted curves correspond to A — —5.05 • 10~ 3 , A = 
and A — —0.012 respectively. Insets in (a) and (b) show 
the real and imaginary parts of permittivity of host material 
respectively. 



region is finite). This discontinuity 20 is accompanied by 
enhanced reflection (i? > 1), and has a physical origin 
similar to the one of enhanced reflectivity reported 
m i6,20 f or g am media excited by evanescent waves in 
total internal reflection geometry. 

To conclude, we developed a universal approach to de- 
termine the sign of refraction index and impedance in 
active and passive media. We have further utilized this 
approach to demonstrate versatile dispersion manage- 
ment, achieving dn p h/dLU = and dZ/du; = regimes in 
nanowire-based NIM system with bandwidth equivalent 
to pico-second optical pulses. The developed technique 
can be readily utilized to determine sign of refraction 
index in different classes of materials and structures, in- 
cluding split-ring- and fish-scale geometries, waveguides- 
based and anisotropy-based NIMs, and can be used to op- 
timize the dispersion of these structures for various pho- 
tonic applications. Furthermore, a combination of sev- 
eral dopants with tailored gain/ absorption spectra can be 
used to engineer a metamaterial having drip^/duj = and 
dZ/doj = in the same frequency range. We note that 
handedness and dispersion of modes in nm-thick meta- 
materials will be strongly affected not only by dispersion 
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of their constituents, but also by overall geometry^!. In 
general, NIM structures with smaller losse o 19 ' 25 will ex- 
hibit weaker dispersion, potentially increasing operating 
bandwidth, and simultaneously reducing the gain values, 
required to achieve efficient dispersion management. In 
particular we expect gain < 100cm -1 to be sufficient in 
nonmagnetic anisotropy-based NIMs^. 
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